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1. Introduction 

A study of the fission barrier heights Bf of nuclei is motivated by the 
importance of this quantity for several physical phenomena. For example, 
the r— process of stellar nucleosynthesis depends (among other quantities 
such as masses and /3-decay rates) on the fission barriers of very neutron- 
rich nuclei 1,2 . The population and survival of hyperdcformcd states at high 
spin also depends on the fission barriers 3 ' 4 . In addition, the physics of fission 
barriers is intimately connected with on-going search for new superheavy 
elements (SHE). The probability for the formation of a SHE in a heavy- 
ion- fusion reaction is directly connected to the height of its fission barrier 5 ; 
the large sensitivity of the cross section a for the synthesis of the fissioning 
nuclei on the barrier height Bf also stresses a need for accurate calculations 
of this value. The survival of the actinides and SHE against spontaneous 
fission depends on the fission barrier which is a measure of the stability of 
a nucleus reflected in the spontaneous fission lifetimes of these nuclei 7 . 

The recent progress in the description of fission barriers within covari- 
ant density functional theory (CDFT 6 ) is briefly reviewed in the current 
manuscript. The comparison with the results of other model calculations is 
presented. Open problems in the study of fission barriers are also discussed. 
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Fig. 1. The difference between experimental and calculated heights of inner fission 
barriers as a function of neutron number N. The results of the calculations are compared 
to estimated fission barrier heights given in the RIPL-3 database 8 , which is used for this 
purpose in the absolute majority of theoretical studies on fission barriers in actinides. The 
results of the calculations within microscopic+macroscopic method ('MM(Dobrowolski)' 
9 and 'MM (Moller)' 10 ), covariant density functional theory ('CDFT' xl ) and density 
functional theory based on the finite range Gogny force ('Gogny DFT' 12 ) are shown. 
Thick dashed lines are used to show the average trend of the deviations between theory 
and experiment as a function of neutron number. The average deviation per barrier 8 [in 
MeV] is defined as S = ^2 i=1 \B\(th)—B\(exp)\/N, where N is the number of the barriers 
with known experimental heights, and B^(th) (B^exp)) are calculated (experimental) 
heights of the barriers. Long-dashed lines represent the trend of the deviations between 
theory and experiment as a function of neutron number. They are obtained via linear 
regression based on a least square fit. From Ref. 13 . 



2. Current status of the investigation of fission barriers in 
actinides 

The progress in the development of computer codes and the availability 
of powerful computers has allowed to study in a systematic way the ef- 
fects of triaxiality on the fission barriers leading to their realistic descrip- 
tion. Within the CDFT framework, the inner fission barriers with triax- 
iality included have been studied for the first time in Ref. 11 using tri- 
axial RMF+BCS approach and the NL3* parametrization 14 . Two years 
later, similar studies of fission barriers in actinides have been performed in 
Refs. 15 < 16 using the RMF+BCS framework with the PC-PK1 and DD-PC1 
paramctrizations of the RMF Lagrangian. The accuracy of the description 
of the heights of inner fission barriers in these calculations is comparable 
with the one obtained in Ref 11 . The calculations of Refs. 15 ' 16 also include 
the results for outer fission barriers where the effects of octupole deforma- 
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tion (and triaxiality [Ref. 15 ]) are taken into account. They agree reasonably 
well with experimental data. 

As compared with axially symmetric calculations, the inclusion of tri- 
axiality has drastically improved the accuracy of the description of inner 
fission barriers in all model calculations (see Ref. 13 and references quoted 
therein). Fig. l a shows that the state-of-the-art calculations within different 
theoretical frameworks (including CDFT) are characterized by comparable 
accuracy (the <5-values) of the description of inner fission barriers. Note that 
this comparison covers only the results of systematic triaxial calculations of 
even-even Th, U, Pu, Cm and Cf nuclei. Recent calculations within Skyrme 
DFT are also characterized by similar accuracy of the description of fission 
barriers 17 ' 18 . 

Minor differences between the approaches in the obtained average devi- 
ations per barrier (Fig. 1) are not important considering the considerable 
uncertainties in the extraction of inner fission barrier heights from experi- 
mental data (see discussion in Ref. 19 ). However, the similarity of the av- 
erage trends of these deviations (shown by thick dashed lines in Fig. 1) as 
a function of neutron and proton numbers is more important considering 
the differences in underlying mean fields and in the treatment of pairing 
correlations. At present, no clear explanation for these trends is obtained 

13,19 

The investigations of Refs. n > 15 > 16 clearly indicate that good description 
of fission barriers in actinides in the CDFT framework is obtained with 
the parametrizations fitted with no information on fission barriers. On the 
contrary, successful description of fission barriers in actinides within non- 
relativistic DFT is based on the parametrizations which explicitly use either 
fission barrier heights (SkM* in Skyrme DFT and D1S in Gogny DFT) or 
fission isomer excitation energies (UNEDF1 in Skyrme DFT); the latter 
being strongly correlated with inner fission barrier height. 

3. Fission barriers in superheavy nuclei. 

Fig. 2 shows how the models which have been bench-marked in a systematic 
way in the actinides (see Fig. 1) extrapolate to the region of superheavy 
nuclei. Note that the results of the MM calculations of Ref. 20 labeled as 
'MM(Kowal)' are not shown on Fig. 1. However, they describe inner fission 
barriers of actinides very accurately. One can see that the model predictions 
vary wildly; the difference in inner fission barrier heights between different 



a Similar figure but as a function of proton number Z is shown in Ref. 
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Fig. 2. Inner fission barrier heights Bf as a function of neutron number N. The results 
of the MM calculations are taken from Ref. 10 (labeled as 'MM (Moller)') and Ref. 20 
(labeled as 'MM (Kowal)'). The position of inner fission barrier saddle in deformation 
space varies as a function of particle number. The labeling of Ref. 13 is used is order 
to indicate whether the saddle is axial (labeled as 'Axial'), has small (7 ~ 10°, labeled 
as 'Ax-Tr') or large (7 ~ 25°, labeled as 'Tr-A') 7-deformations in the RMF+BCS 
calculations. 

models reaches 6 MeV in some nuclei. This is despite the fact that these 
models describe the inner fission barriers in actinides with comparable level 
of accuracy. The more surprising fact is that the prediction of two MM mod- 
els differ so substantially; in reality the 'MM (Kowal)' model predictions 
are closer to the CDFT ones than to the 'MM (Moller)' predictions. 

4. Some open questions 

The results in the SHE region raise a number of questions which require 
further investigation. Among these, for example, are 

• A. How the pairing evolves with deformation? How well the pairing 
strength is defined for SHE? 

• B. How the accuracy of the description of the energies of the single- 
particle states and its evolution with deformation affects the results 
of the calculations? 

I will concentrate on question A. Table IV in Ref. 13 shows the large 
variety of pairing prescriptions used in the calculations of fission barriers 
during last decade. The question definitely emerges how well they describe 
the evolution of pairing with deformation and extrapolate towards SHE. 
These features are important for a quantitative description of fission barrier 
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Fig. 3. Experimental and calculated neutron three-point indicators A^(JV) as a func- 
tion of neutron number N. The results of the CRHB and CRHB+LN calculations with 
the NL1 parametrization are shown. 



heights which according to Ref. 23 sensitively depend on pairing properties. 
Although some attempts were made in 70ies to extract the information 
on pairing properties at fission saddles of actinides 24 , they did not lead 
to reliable estimates. It is even more difficult to get information on pair- 
ing properties of SHE due to the lack of relevant experimental data. As 
a consequence, as accurate as possible description of pairing properties in 
actinides and light SHE should be considered as a necessary condition for 
further extrapolation to heavier SHE. 

In order to address these questions, a comprehensive study of pairing 
and rotational properties of actinides is performed in Ref. 25 in the cranked 
relativistic Hartree-Bogoliubov (CRHB) framework. In this approach, the 
Brink-Booker part of phenomenological non-relativistic Gogny-type finite 
range interaction 

V pp (l,2) = f e~ [{ri ~ r2), ^ ? {Wi + B t P a - H,^ - NhP 17 P T ) (1) 

i=l,2 

is used in pairing channel. In addition, an approximate particle number pro- 
jection by means of Lipkin-Nogami (LN) method is employed. The analysis 
indicates that the pairing strength has to be decreased as compared with 
original D1S Gogny force; the scaling factor / is equal to 0.899 and 0.9147 
for the NL3* and NL1 parametrizations, respectively. Only with these scal- 
ing factors, the moments of inertia and odd-even mass staggerings (the A^ 3 ) 
indicators) can be well described in the normal-deformed (ND) minimum 
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Fig. 4. Calculated and experimental moments of inertia at low spin. 



(Figs. 3 and 4). 



Table 1. Experimental and theoretical charge quadrupole moments Q of SD fission iso- 
mers. The results of the CRHB+LN calculations with the NL1 and NL3* parametriza- 
tions are presented. Experimental data for the U and Pu isotopes are taken from Ref. 
21 , while the one for 242 Am from Ref. 22 . 
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Fission isomers provide only available tool to estimate the evolution 
of pairing with deformation in actinides. Such an estimate is available 
only through the study of rotational properties of 236 > 238 TJ and 240 Pu 
nuclei; these are only nuclei for which supcrdcformcd (SD) rotational 
bands were experimentally measured. The calculated charge quadrupole 
moments Q are compared with available experimental data in Table 1. 
The CRHB+LN(NL3*) results of the calculations come reasonably close to 
experiment. The CRHB+LN(NL1) results are also not far away from ex- 
perimental data but they substantially overestimate experimental Q value 
in 238 U. 

The experimental kinematic moments of inertia are best described by 
the NL3* parametrization; the deviation from experiment does not exceed 
3.4% (Fig. 5). The fact that the moments of inertia of rotational structures 
in two different minima (ND and SD) are accurately described with the 
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same pairing strength strongly suggests that the evolution of pairing cor- 
relations with deformation is properly described in the CRHB+LN(NL3*) 
framework by the Gogny D1S pairing. 
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Fig. 5. Experimental and calculated kinematic (J^ 1 ') and dynamic (J^ 2 ') moments of 
inertia of SD rotational bands in 236 > 238 TJ and 240 Pu. The notation of the lines and 
symbols is given in the figure. 

However, this is not always the case since the CRHB+LN(NL1) cal- 
culations substantially overestimate the experimental moments of inertia 
in the SD minimum [for example, by 11.3% in 240 Pu] (Fig. 5) while they 
reproduce the low-spin moments of inertia in the ND minimum with the 
same level of accuracy as the CRHB+LN(NL3*) calculations (see Fig. 4). 
It turns out that reasonable description of the moments of inertia at SD 
can be achieved in the CRHB+LN(NL1) calculations only if scaling factor 
/ « 1.0 is used. This clearly indicates that even the pairing force carefully 
fitted to experimental data at ND does not guarantee accurate description 
of pairing at SD (and as a consequence also at fission saddle). 

The origin of such behavior is not completely clear but the difference be- 
tween the CRHB+LN(NL3*) and CRHB+LN(NL1*) results for jW at SD 
may also partially originate from the differences in the single-particle struc- 
tures at superdeformation obtained with the NL3* and NL1 parametriza- 
tions. The Q values obtained in the CRHB+LN(NL1) calculations are al- 
ways higher than the ones for CRHB+LN(NL3*), which also may be a 
reason why the CRHB-f-LN(NLl) calculations systematically overestimate 
kinematic moments of inertia at SD. 

5. Conclusions 

The review of the recent applications of covariant density functional theory 
to the description of fission barriers in actinides and superheavy nuclei is 
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presented. Inner fission barriers of actinides are described with compara- 
ble accuracy in the MM and DFT (including CDFT) calculations. How- 
ever, the predictions for fission barriers of SHE within different frameworks 
vary drastically indicating the need for better parametrizations of different 
channels of the models. In particular, a better pairing force and its better 
paramctrization is needed. An effort to improve the Gogny pairing force 
in the RHB framework, the first results of which are presented here, is 
currently in progress. 
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